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The problem of parameter estimation on hybrid- wing-body type aircraft is complicated 
by the fact that many design candidates for such aircraft involve a large number of aero- 
dynamic control effectors that act in coplanar motion. This fact adds to the complexity 
already present in the parameter estimation problem for any aircraft with a closed-loop 
control system. Decorrelation of system inputs must be performed in order to ascertain 
individual surface derivatives with any sort of mathematical confidence. Non-standard 
control surface configurations, such as clamshell surfaces and drag-rudder modes, further 
complicate the modeling task. In this paper, asymmetric, single-surface maneuvers are used 
to excite multiple axes of aircraft motion simultaneously. Time history reconstructions of 
the moment coefficients computed by the solved regression models are then compared to 
each other in order to assess relative model accuracy. The reduced flight-test time required 
for inner surface parameter estimation using multi-axis methods was found to come at the 
cost of slightly reduced accuracy and statistical confidence for linear regression methods. 
Since the multi-axis maneuvers captured parameter estimates similar to both longitudinal 
and lateral-directional maneuvers combined, the number of test points required for the 
inner, aileron- like surfaces could in theory have been reduced by 50%. While trends were 
similar, however, individual parameters as estimated by a multi-axis model were typically 
different by an average absolute difference of roughly 15-20%, with decreased statistical sig- 
nificance, than those estimated by a single-axis model. The multi-axis model exhibited an 
increase in overall fit error of roughly 1-5% for the linear regression estimates with respect 
to the single-axis model, when applied to flight data designed for each, respectively. 


Nomenclature 


a angle of attack 

q dynamic pressure 

/ 3 angle of sideslip 

5a antisymmetric (roll) paired surface deflection 

5e symmetric (pitch) paired surface deflection 

5r surface deflection of winglet rudders, positive for postive yaw 

5s j deflection of surface j, positive for trailing edge down 

0 vector of model parameters 

b reference span 

C a coefficient of axial force, 

Ci coefficient of rolling moment, ||| 

Cm coefficient of pitching moment, 

C/v coefficient of normal force, ^ 

C n coefficient of yawing moment, ||| 
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Cy coefficient of side force, 

C FSaj derivative of force or moment coefficient [F] with respect to antisymmetric surface pair [j] deflection 

Cpbej derivative of force or moment coefficient [F] with respect to symmetric surface pair [j] deflection 

Cpssj derivative of force or moment coefficient [F] with respect to asymmetric surface [j] deflection 

Cmoi derivative of coefficient of pitching moment with respect to angle of attack 

Cno. derivative of coefficient of normal force with respect to angle of attack 

C na derivative of coefficient of yawing moment with respect to angle of sideslip 

Cyp derivative of coefficient of side force with respect to angle of sideslip 

g acceleration due to gravity 

I identity matrix 

I xx roll moment of inertia 

I xz XZ-plane cross-axis moment of inertia 

I yy pitch moment of inertia 

I zz yaw moment of inertia 

1 when applied to a surface deflection, a subscript to indicate the left-side surface 

M Fisher information matrix 

m aircraft mass 

p roll rate 

r yaw rate; or, when applied to a surface deflection, a subscript to indicate the right-side surface 

r 2 statistical coefficient of determination 

S reference area 

s regression model fit error 

s 2 parameter variance 

V magnitude of flight velocity 

X regressor matrix 

2 output vector for the linear regression problem 
AFRL Air Force Research Laboratory 

DFRC Dryden Flight Research Center 
HWB Hybrid-Wing-Body (class of aircraft) 

NASA National Aeronautics and Space Administration 

I. Introduction 

The X-48B Blended Wing Body is an 8.5% dynamically-scaled, hybrid-wing-body (HWB) aircraft built 
by Cranfield Aerospace Ltd (United Kingdom); it was the result of a joint partnership between the National 
Aeronautics and Space Administration (NASA), the Air Force Research Laboratory (AFRL), and The Boeing 
Company (Chicago, Illinois). The X-48B and its planned successors are representative of possible future, 
highly efficient, HWB transport designs that involve integrated configurations using a large number of control 
effectors. The Environmentally Responsible Aviation project within the NASA Aeronautics Research Mission 
Directorate intends to use the conclusions drawn from flight research on a series of configurations for the 
X-48B to study the development of new, manned, highly-efficient, HWB transport aircraft. A photograph 
of the X-48B in flight at the NASA Dryden Flight Research Center (DFRC) (Edwards, California) is shown 
in Fig. 1. 

Research is being conducted at DFRC to study the effectiveness of various parameter estimation methods 
on HWB aircraft. The problem of parameter estimation on HWB aircraft is complicated by the fact that 
many design candidates for such aircraft involve a large number of aerodynamic control effectors that act 
in coplanar motion. This adds to the complexity already present in the parameter estimation problem for 
any aircraft with a closed-loop control system. Decorrelation of system inputs must be performed in order 
to ascertain individual surface derivatives with any sort of mathematical confidence. Non-standard control 
surface configurations, such as clamshell surfaces and drag-rudder modes, further complicate the modeling 
task. 

A previous paper 1 that studied the efficacy of single-axis, doublet-based parameter estimation maneuvers 
in the longitudinal axis was published by two of the present authors in early 2010. A paper 2 on the lateral- 
directional results using these same methods was published in 2011. These two papers establish a set of 
results using single-axis methods that were used as the baseline for the present study. 
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Figure 1. The Boeing X-48B Blended Wing Body, shown in flight near the NASA Dryden Flight Research 
Center, Edwards, California. 


The present paper focuses on indentifiability considerations of performing simultaneous parameter esti- 
mation for all three axes of motion using a single-surface maneuver. Such maneuvers will be referred to as 
“multi-axis” maneuvers. Multi-axis maneuvers carry the benefit of decreased flight-test time required, since 
multiple parameters may be solved for using the same maneuver. The tradeoff between reduced flight-test 
time and model accuracy is assessed in this study for the inner, elevon-like surfaces of the X-48B Blended 
Wing Body. 

The nature of the tradeoff between flight-test time and parameter estimation quality may be of interest 
to flight programs for which the amount of total testing is constrained below technically desirable levels by 
financial or operational considerations. 

II. Aircraft Description 

The X-48B aircraft incorporates a unique configuration and outer mold line. Instruments relevant to 
parameter estimation include dual airdata probes to measure airspeed, angle of attack, and angle of sideslip. 
Additionally, the aircraft is equipped with an Inertial Measurement Unit (IMU) and Global Positioning 
System (GPS) that provides linear acceleration, angular rotation rates, Euler angles, and position. Twenty 
actuating aerodynamic surfaces, eighteen of which are coplanar, are used to provide aircraft control. 

The X-48B aircraft can be configured with the leading edge slats extended or retracted; however, they 
cannot be adjusted in flight. The center of gravity can be adjusted on the ground between forward and 
aft configurations. Allocation of the control surfaces is depicted in Fig. 2 with surface pairs numbered for 
reference. The inner surfaces (1 through 5) are elevons. Surfaces 6 and 7 are split ailerons, or “clamshell” 
surfaces: the top and bottom surface can be moved together to produce roll moments or they can be split 
to produce a yaw moment through differential drag. Rudders are incorporated into the winglets to provide 
additional yaw control and stability. 

Control surface positions are inferred from the measured actuator position and are not measured directly. 
The control surface actuation on the X-48B aircraft consists of an electro-mechanical servo that moves the 
control surface through a linkage. Position measurement is taken at the output shaft of the servo; thus, 
differences between the surface position and actuator position may be due to linkage bending or gear slop. 
No corrections were made to the control surface data because data or models necessary for corrections were 
not available. 

While these factors complicate any attempt to make absolute quantification of the individual parameters, 
much work can be done in studying the relative effectiveness of parameter estimation techniques and the 
associated model structure. Present experimental work by NASA is ongoing to attempt to quantify the 
uncertainty in the control surface position, as well as the propagated effect into the estimated parameters. 
This paper applies the models suggested by stepwise regression to simulation and flight data for the X-48B 
aircraft, analyzes trends in the Cramer- Rao lower bounds, assesses the fit of state- variable time history re- 
constructions from the solved parameters, and verifies the stepwise regression model selection by comparison 
to possible alternatives. 
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Figure 2. Control surface numbering and axis allocation for the X-48B Blended Wing Body. 


III. Method 

Parameter estimation is a subset of the broader field of system identification, wherein the basic task of 
the engineer is to determine the nature of a system under study through observation and analysis of the 
outputs generated by a controlled set of inputs . 3 Parameter estimation assumes that the system in question 
may be modeled as a parametrized set of equations, of which the coefficients, or parameters, are the objective 
of the analysis. 

The generally accepted standard method of performing maneuvers for parameter estimation are the 
traditional doublet or triplet inputs. The pilot inputs a simple square wave command of controlled magnitude 
in a particular axis, and the output dynamics that result from this input are analyzed for a mathematical 
relationship. In the case of multiple surfaces affecting the same axis, in theory, each surface and associated 
surface (its opposite pair on the other wing, as well as adjacent surfaces that may provide interference effects) 
must be tested in all possible combinations for comprehensive model validation. On aircraft with a high 
number of surfaces, this can be quite time-consuming. 

A. Linear Regression for Parameter Estimation 

Linear regression parameter estimation is a technique in which the coefficients of an assumed linear relation- 
ship between known inputs and observed outputs are estimated using least-squared fits. More information 
on linear regression parameter estimation techniques can be found in Klein and Morelli . 4 Linear regression 
is often compared to output-error techniques, which actually integrate the equations of motion and compare 
the state outputs to measured values in order to estimate the parameters. In this sense, linear regression 
is mathematically simpler and therefore quicker in terms of computational time and complexity. The full 
derivation of such techniques is beyond the scope of this paper. More information on output-error and max- 
imum likelihood techniques can be found in many available references on the subject . 3,5-8 Similarly, there 
are many references 9-15 that describe prior work to decorrelate control surfaces, including substantial work 
by Morelli and his colleagues in the field of optimal input design. 

A tool often used to quantify the relative statistical confidence of a parameter estimate is the Cramer- 
Rao lower bound. The Cramer-Rao lower bound represents the lowest magnitude limit for the variance of 
an estimator with a given bias 16 . Un-modeled dynamics can make the true value of the variance of the 
estimator much higher. In the simplest case where the variance is assumed to be unbiased and have a normal 
distribution, the Cramer-Rao bound becomes simply the inverse of the Fisher information matrix, M(^), 
which is a metric for measuring the amount of usable information content in a set of data. Choosing input 
design methods and flight-test techniques that lower the Cramer-Rao bounds is an effective approach to 
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choosing inputs that maximize the usable information content of the flight data. 

The derivation of the Cramer-Rao bound assumes that the residuals consist solely of white noise. In 
practice, this is not the case. As a result, traditionally-computed bounds can be inaccurate. The Cramer- 
Rao bounds presented in this report have been adjusted to account for frequency content in the residuals, 
using a technique from Klein and Morelli. 4 

B. Data Sources 

Time histories of the aircraft angular rates, air data, control surface position, and other pertinent information 
for performing parameter estimation problems were obtained from two sources: the nonlinear simulation 
provided by The Boeing Company (Chicago, Illinois), and recorded data from the Phase 1 and Phase 1.5 
series of flight tests of the X-48B aircraft at NASA DFRC. For both simulation and flight, force and moment 
coefficients were constructed from observed air data and flight dynamics using standard aircraft equations 
of motion. 

A. Maneuver Description 

For the standard, single-axis maneuvers, a series of symmetric (equal magnitude in the same direction, 
inducing pure pitch) or antisymmetric doublets (equal magnitude in opposite directions, inducing pure roll) 
was performed in both simulation and flight. For the simulation studies, the surface motions could be 
performed serially in a single combined maneuver termed a supermanuever, on which regression could be 
performed. An example is shown in Figure 3. 



Figure 3. Multi-surface supermaneuver for temporal separation. 


It is worthwhile to note that the use of supermaneuvers means that every surface excitation in a sequence 
is present in the dataset for all parameter estimate regression solutions, regardless of which parameters are 
the focus of the particular model being solved. For more details on the single-axis parameter estimation 
method and results for the X-48B aircraft, see Taylor and Ratnayake 1 for the longitudinal results, and 
Ratnayake, Waggoner, and Taylor 2 for the lateral-directional results. 

For multi-axis / single-surface maneuvers, similar supermanuevers were constructed for the inner five 
surfaces, but using only a surface on one wing instead of using paired symmetric or antisymmetric maneuvers. 
This asymmetry induces motion in multiple axes simultaneously. Parameter estimation for all surfaces under 
consideration was then performed on the single combined supermaneuver. 

B. Flight- Testing 

Flights 65, 66, and 67 from Phase 1 of the X-48B flight-test program included standard, single-axis, doublet 
lateral-directional parameter estimation test points for the present analysis, at angles of attack of 6 deg, 8 
deg, and 10 deg, respectively. Flights 85 and 89 from Phase 1.5 of the program included single-surface / 
multi-axis maneuvers at the same angles of attack. The range of altitudes and corresponding trim speeds 
for these angles of attack do not vary enough for a second variable to be considered in defining the flight 
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condition. Individual surface pair maneuvers were available for analysis. Each maneuver was performed 
three times in flight in order to reduce error and quantify the variance of the dataset. 

While smoothing was applied to the transition points to eliminate any discontinuities, two additional 
factors ameliorated the effect of splicing on the parameter estimation results. First, linear regression, as an 
equation-error technique, does not rely on integrating the equations of motion as in output-error techniques. 
This means that there is no need to ensure that integrators are properly reset and that data is exactly 
aligned at the splicing points for the parameter estimation (though it is required later for time-history 
reconstruction). Second, the range of flight conditions and possible trim states for the X-48B aircraft is 
limited enough that significant disparities in state variables for the same trim angle of attack are unlikely. 

Noise present in the flight data was removed using a third-order, two-way Butterworth filter applied with 
a corner frequency of 5 Hz, because this cutoff appeared to capture the dynamics of interest while excluding 
the bulk of the measurement and physical noise. Corrections were also made to measured air data as well 
as translational accelerations to account for the distance from the aircraft center of gravity to the respective 
measurement points. The airdata is thus corrected for upwash and sidewash induced by the rotation rates 
of the aircraft. 

C. Nonlinear Simulation 

The nonlinear simulation of the X-48B aircraft was designed by The Boeing Company and is implemented 
in Simulink@ (The MathWorks, Natick, Massachusetts). The version of the simulation used for this paper 
was 4.3.1, using Vehicle Management System (VMS) version 4.3 and aerodynamic model 20100527. Note 
that this is a more recent aerodynamic model version than used in previous work that was performed on 
Phase 1 data only . 1,2 

Though flight is the only aeronautical reality, the advantage of conducting simulation studies is that 
various real-world effects can be controlled or are often simply not modeled, and the isolated effects of 
various external factors on the results can be explored. The simulation results shown in this paper were for 
supermaneuvers performed with the flight control laws forced into an open-loop mode. Further, the normal 
control surface allocator was bypassed in order to insert specific input combinations that were not available 
in the normal control laws. Applying the parameter estimation method on the simulation data also allows 
for consistency checks on the method. 

The simulation maneuvers were initially planned for the same angles of attack as the flight data. It 
became necessary, however, to move the simulation points to take place at trim angles of attack on the half- 
angle (for example, 6.5 deg, 7.5 deg, et cetera) in order to avoid breakpoints in the aerodynamic model, which 
were presumably artifacts from the wind-tunnel test procedure. These breakpoints made finite difference 
approximations in their vicinity difficult, which affects related research that would use the same data. 

C. Effector Definitions 

The configuration of the X-48B aircraft allows for several possible definitions of the control effector regressor 
functions due to the split nature of the clamshell surfaces. Because the clamshell surfaces can move in unison 
or oppositely, the same clamshell upper and lower surface can behave like a traditional aileron (primarily 
affecting roll), or split open in a drag-rudder yaw mode. 

Deflection of any individual surface will be denoted as Ss; for example, the deflection of the inner surfaces 
will be represented by through Ss 5 , where the number corresponds to the surface number shown in Fig. 
2 . When treated individually, the inboard clamshell upper and lower surfaces will be denoted as Ssq u and 
Ssqi , respectively. The outboard clamshell upper and lower surfaces will be similarly denoted Ss^ u and Ssn. 
The winglet rudders are located on the wingtips and behave as normal rudders (affecting the directional 
axis only). The deflection of the winglet rudders is denoted Sr. When an inner surface is coupled into 
an anti-symmetric pair with its counterpart on the opposite wing of the aircraft, the aileron-like combined 
deflection can be described by a single abstracted effector, Sa , as shown in Eq. ( 1 ) . A similar definition for 
symmetric, elevator-like motion is described by Eq. (2). 

— $Sj,left right ( 1 ) 

Sej = 8sjj e ft T Ssj, r ighi (2) 
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D. Model Structure 

The multi-axis parameter estimation problem seeks to find the linear parameters that, when multiplied by 
their respective regressor functions, compute accurate values for axial force (Ca), normal force (Cjy), side 
force (CV), rolling moment (C|), pitching moment (C m ), and yawing moment ( C n ). The assumed relationship 
between the regressor functions and the output force and moment coefficients is the model structure. The 
primary metrics for comparing the quality of the parameter estimates made by traditional methods versus 
multi- axis methods are the variance of the estimates, the Cramer- Rao lower bounds, and the time history 
reconstructions. 

The multi-axis maneuvers are analyzed by solving the parameter estimation problem for both longitudi- 
nal and lateral-directional parameters simultaneously. The results of these analyses are then compared to 
parameter estimates obtained from single- axis maneuvers, which are analyzed by solving the equations for 
that axis only. The comparison then is between the single-axis (traditional) maneuvers which isolate motion 
as much as possible into ideally pure roll, pitch, or yaw, and multi-axis maneuvers which attempt to excite 
all axes at the same time. 

For linear regression analysis, the general form is shown in Eq. (3). 


z = X0, (3) 

In Eq. (3), 2 is the column array of force and moment coefficients, X is the matrix of regressor functions 
(with individual model equations arranged in rows), and 0 is the column array of parameters to be estimated. 
The force and moment coefficients are calculated from known or measured aircraft dynamics such as linear 
and rotational acceleration, instantaneous center of gravity, and moments of inertia. The regressor functions 
are known or constructed from aircraft state measurements, such as the angles of attack and sideslip, surface 
positions, and the instantaneous dynamic rates. Solving for the parameters yields Eq. (4): 


0 = (. X T X)~ 1 X T z (4) 

For single-axis maneuvers, these equation sets are solved for in-axis derivatives only, in response to pure 
axis inputs. For example, a left surface and its corresponding right surface moving anti-symmetrically would 
produce (ideally) a pure roll input, and a response from the aircraft that should be dominated by roll-axis 
dynamics. An example equation for coefficient of rolling moment would be Eq. (5): 

Cl = Cio + Cl p • /? + Cip • + Clr • + [Cl ^surfaces] (5) 

The first four terms of Eq. 5 are quite straightforward; the rolling moment should depend on some bare- 
airframe damping term (the regressor function for which is simply 1), the angle of sideslip, and the roll and 
yaw rates of the aircraft. These aerodynamic stability and damping terms will be combined and referred to 
as Ci, aero (and like manner for the directional coefficients) in the remainder of the paper; however, it should 
be noted that the four described components of Ci, aero are solved for individually as their own regressor 
functions. An analogous example can be inferred for symmetric maneuvers acting purely on the pitch axis, 
or symmetric rudder maneuvers acting purely on the yaw axis. 

In the case of single-surface, multi-axis excitations, the inputs are neither symmetric nor antisymmetric. 
Instead, a single surface on one side of the aircraft is subjected to a doublet input, and the asymmetry 
results in an excitation in the pitch, roll, and yaw axes simultaneously. Thus, the rolling moment coefficient 
equation becomes Eq. (6): 

Cl, aero = Qo + Ci a • O' + Cip * /3 + Ci p • + Ci q • + Ci r • (6) 

In Eq. 6, longitudinal regressors for angle of attack and pitch rate are added to capture multi-axis 
dynamics. The corresponding C/ ?sur / aces term would then consist of the the individual surface motions, as 
described by Eq. 7. 


Ci ,surfa 


. Cl$Sn ■ djSj 


(7) 


The parameters are solved for using analogous equations for the remaining force and moment coefficients 
at each time frame of the data. The measured input data can be passed through the resulting equations (using 
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the estimated parameters) to arrive at the equations of motion for the aircraft, which are then integrated to 
arrive at time history reconstructions of key state variables, such as angle of attack. These reconstructions 
can be compared to the measured, real-time histories to assess the accuracy of the parameter estimates. See 
Klein and Morelli 4 or Maine and Iliff 7 for more detailed explanations of reconstructing the time histories 
using aircraft equations of motion, as well as the formulation of the force and moment coefficients used in 
this paper. 

E. Statistical Analyses 

Statistical analyses were performed on representative samples of flight data to determine the relative expected 
performance of the single-axis versus the multi-axis model. For the overall model equations, the tools used 
were the model coefficient of determination (r 2 ) and the fit error (s), defined in Eqs. (8) and (9), respectively. 
The average results for single- and multi-axis flight data taken at an angle of attack of 6 deg are shown in 
Table 1. 


2 _ & T X T z - Nz 2 
z T z — Nz 2 


(8) 


s = 


EQi - in ) 2 

N — p 


(9) 


Model 

Structure 

r 2 

s 

Cl 

Single-axis 

92.10% 

27.80% 


Multi-axis 

90.66% 

30.36% 

c m 

Single-axis 

93.83% 

24.77% 


Multi-axis 

92.69% 

26.97% 

C n 

Single-axis 

90.44% 

30.79% 


Multi-axis 

90.42% 

31.50% 


Table 1. Average coefficient of determination (r 2 ) and fit error (s) for single-axis and multi-axis model struc- 
tures, performed on three repetitions of representative flight data at 6 deg trim angle of attack. 


It should be noted that the results in Table 1 are for each model as applied to a supermaneuver constructed 
from sub-maneuvers that were designed for the model in question. In other words, the lateral-directional 
single- axis model was applied to a supermaneuver of anti-symmetric doublets, the longitudinal single-axis 
model was applied to a supermaneuver of symmetric doublets, and the multi- axis model was applied to a 
supermaneuver of asymmetric, single-surface doublets. 

The coefficient of determination and the fit error are metrics that apply to the model equation as a 
whole. In addition, the Fo-statistic may be used to assess the relative statistical significance of of individual 
parameters in the regression fit. The Fo-statistic is described in Eq. (10), in which s 2 in this case is the 
parameter variance [not explicitly related to the fit error in Eq. (9)]. 


F 0 = 



( 10 ) 


Tables 2 and 3 show the F-ratio for rolling and pitching moment derivatives, respectively; they were 
calculated using the same method-specific samples of flight data as in the analyses of coefficient of determi- 
nation and the fit error. In these tables, a higher value of the F-ratio connotes a higher level of significance 
to the regression fit. 

For both the roll and the pitch axes, the multi- axis model is shown to exhibit F-ratios for the surface 
derivatives that are of the same order of magnitude as the single-axis model, though notably less in value. 
Overall, the significance values appear high; it should be noted that the manuevers in the flight data used 
for the analysis were designed specifically for these models. 

All three statistical analyses suggest a priori that multi-axis models are likely to perform slightly worse 
than single-axis models with respect to model accuracy. 
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Cl 

Sa 

<5s 

Sfc 2 

3.51 • 10 4 

1.25 • 10 4 

Sfc 3 

4.58 • 10 4 

1.69 • 10 4 

Sfc 4 

3.24 • 10 4 

1.63 • 10 4 

Sfc 5 

4.89 • 10 4 

1.74 • 10 4 


Table 2. Partial F-ratio for derivative of coefficient of rolling moment with respect to single-axis (Sa) and 
multi-axis (5s) surface deflection. 


Cm 

Se 

Ss 

Sfc 2 

2.19- 10 4 

1.85 • 10 4 

Sfc 3 

2.28 • 10 4 

1.83 • 10 4 

Sfc 4 

2.39 • 10 4 

1.65 • 10 4 

Sfc 5 

2.41 • 10 4 

1.63 • 10 4 


Table 3. Partial F-ratio for derivative of coefficient of pitching moment with respect to single-axis (Se) and 
multi-axis (Ss) surface deflection. 


IV. Results and Discussion 

This section presents results of the comparison between the single- axis and multi- axis methods. Due to 
the proprietary nature of the performance data of the X-48B aircraft, all plots in this section are provided 
without quantification on the ordinate axis. The flight-data results are provided with error bars, which 
represent the Cramer-Rao lower bounds for the respective parameter estimate data points. 

Many of the provided results depict trends for simulation and flight that do not agree with each other in 
slope or other characteristics. While the simulation and flight results will be loosely compared as a sanity 
check to ensure that the parameter estimation results for flight do not wildly deviate from expected values, 
it is not the objective of this paper to assess the accuracy of the Boeing nonlinear simulation or aerodynamic 
model with respect to flight. The objective is rather to compare the accuracy of the single-axis and multi-axis 
model on multiple HWB aircraft datasets. 

The best measure of the effectiveness of a model lies in its ability to more accurately reconstruct observed 
state- variable time histories. In doing so, the best model should also perform with lower Cramer-Rao bounds 
than do other models for the same dataset. The models are here applied to two different datasets for the 
X-48B aircraft: simulation and flight. Any disparity between the two sets of results suggests that the 
aerodynamic model inadequately captures certain aerodynamics of the flight vehicle; it does not affect the 
comparison between the parameter estimation models. 

Results are not shown for the innermost surface pair, surfaces 1, since the X-48B flight control laws do 
not permit a lateral roll input using this surface pair; they are designed for pitch control only. 

A. Aerodynamic Derivatives 

Representative aerodynamic derivatives for force and moment in the longitudinal and lateral-directional axes 
are shown in Figs. 4 and 5. 

The multi-axis methods deliver estimates of aerodynamic derivatives that follow similar trends to the 
single-axis estimates. The multi-axis method, however, is shown to exhibit greater spread (corresponding 
to reduced consistency) and wider Cramer-Rao lower bounds (corresponding to increased uncertainty) than 
the single-axis method. Additionally, the average absolute percent difference in the individual parameter 
estimates can to be quite high, as shown in Table 4. 

Cy (3 in particular appears to be estimated poorly. It should be noted, however, that in some cases the 
magnitudes of the derivatives are quite small. Thus, a small-magnitude difference in the estimates can lead 
to relatively high percentage differences. As another perspective on the error, for Cyp, the relative difference 
between methods is approximately 6 times the span of the Cramer-Rao lower bound. The weakness in the 
modeling of Cyp by the multi-axis method is presumed to be related to the unique geometry of the X-48B 
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• Sim SA 
□ FltMA 
A FltSA 


(a) Coefficient of normal force with respect to angle of at- (b) Coefficient of pitching moment with respect to angle of 
tack. attack. 


Figure 4. Representative longitudinal aerodynamic derivatives for single-axis and multi-axis flight data. 



• Sim SA 
□ Fit MA 
A FltSA 



• Sim SA 
□ FltMA 
A FltSA 


(a) Coefficient of side force with respect to angle of sideslip, (b) Coefficient of yawing moment with respect to angle of 

sideslip. 


Figure 5. Representative lateral-directional aerodynamic derivatives for single-axis and multi-axis flight data. 


Derivative 

Diff. 

c Na 

16.78% 

Cma 

18.82% 

Gyp 

82.98% 

C n /3 

13.51% 

c w 

48.45% 


Table 4. Average absolute percent difference between multi-axis and single-axis estimates of selected aerody- 
namic derivatives. 
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aircraft, which is relatively flat in the XY plane; hence it has directional profile which does not change 
significantly with small changes in sideslip angle. 


B. Surface Derivatives 

Example parameter estimates for the pitch and roll effectiveness of the inner surfaces is shown in Figs. 6 
through 9. In these figures, the parameters that were estimated using multi-axis maneuvers (designated 
“MA”) were solved for using regression equations that included the full set of longitudinal and lateral- 
directional terms. The single-axis maneuvers (designated “SA”) were solved for using regression equations 
that isolated each axis, and used data sets from separate single-axis maneuvers. 



• Sim SA 
□ Fit MA 
A FltSA 



• Sim SA 
O' FltMA 
A FltSA 


(a) Roll effectiveness. 


(b) Pitch effectiveness. 


Figure 6. Estimates of surface 2 pitch and roll moment effectiveness. 



• Sim SA 
□ Fit MA 
A FltSA 



• Sim SA 
□ FltMA 
A FltSA 


(a) Roll effectiveness. 


(b) Pitch effectiveness. 


Figure 7. Estimates of surface 3 pitch and roll moment effectiveness. 


Similar to the aerodynamic derivatives, the multi-axis methods deliver estimates of aerodynamic deriva- 
tives that follow similar trends to the single-axis estimates. The average absolute percent difference in the 
individual parameter estimates are shown in Table 5. 


Derivative 

Diff., Sfc 2 

Diff., Sfc 3 

Diff., Sfc 4 

Diff., Sfc 5 

Cl8s 

6.88% 

16.86% 

25.26% 

30.38% 

CmSs 

13.59% 

14.11% 

7.06% 

5.61% 


Table 5. Average absolute percent difference between multi-axis and single-axis estimates of selected surface 
derivatives. 

The pitch axis parameters are in general tightly clustered within each other’s Cramer-Rao lower bounds 
for both single-axis and multi- axis models, though they appear to have larger Cramer-Rao lower bounds 
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Less Effective p More Effective Less Effective p More Effective 



— Sim SA 
□ Fit MA 
A FltSA 



• Sim SA 
□ FltMA 
A FltSA 


(a) Roll effectiveness. 


(b) Pitch effectiveness. 


Figure 8. Estimates of surface 4 pitch and roll moment effectiveness. 




(a) Roll effectiveness. 


(b) Pitch effectiveness. 


Figure 9. Estimates of surface 5 pitch and roll moment effectiveness. 
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overall. The multi-axis estimates for rolling moment surface derivatives, however, agree less with the more- 
established single-axis method as one moves further outward on the wing. The multi-axis estimates for 
pitching moment surface derivatives in general appear to agree more with the single-axis methods as one 
moves further outward on the wing. 

C. Time History Reconstruction 

Time history reconstruction of state variables and moment coefficients is one of the most direct methods of 
assessing the quality of a parameter estimation model. When driven with the same inputs as simulation or 
flight, the model should be able to accurately reconstruct the observed outputs of the system under study. 
The calculated and observed moment equations for simulation and flight, respectively, were constructed using 
Eqs. (11) through (13). 


Q_Sb . I XZ . I ZZ Iyy I XZ 

-j— c i =p - i~ r + — j Q r - -j—qp 

qSc . I xx — I zz I xz , 2 2\ Ip o 

-j—C m = q+ pr + —{p - r ) - -f-Sl p r 

1 yy 1 yy 1 yy 1 yy 

qSb . I XZ . Iyy I xx I xz Ip 

—C n = r- —p + pq + — qr + -f-Q p q 

J-zz J-zz J-zz J-zz J-zz 


( 11 ) 

( 12 ) 

(13) 


The parameter estimates from each model were used to reconstruct estimates of the above experimental 
values using the parameter equations as in Eqs. 5 and 6. The time histories of these were then compared to 
the calculated (simulation) or observed (flight) time histories. The force and moment coefficients can then 
be integrated using the equations of motion to arrive at the aircraft state variables, such as wind angles and 
rotational rates. 


A. Simulation 

Example time history reconstructions of the rolling, pitching, and yawing moments for simulation data are 
shown in Figs. 10 through 12. These figures represent the same slice from a supermaneuver assembled for 
9.5 deg angle of attack. 


At 09p5 deg Trim Alpha At 09p5 deg Trim Alpha 




Figure 10. Coefficient of rolling moment for single-axis and multi-axis models, reconstructed from parameter 
estimates on simulation data at 9.5 deg trim angle of attack. 


The reconstructions of the coefficients on simulation data verify that the models behave as expected 
in a nearly ideal experimental context: open-loop, noiseless, and perfectly controlled. The single-axis and 
multi-axis models reconstruct the moment coefficients nearly perfectly. 
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At 09p5 deg Trim Alpha 


At 09p5 deg Trim Alpha 




Figure 11. Coefficient of pitching moment for single-axis and multi-axis models, reconstructed from parameter 
estimates on simulation data at 9.5 deg trim angle of attack. 


At 09p5 deg Trim Alpha At 09p5 deg Trim Alpha 




Figure 12. Coefficient of yawing moment for single-axis and multi-axis models, reconstructed from parameter 
estimates on simulation data at 9.5 deg trim angle of attack. 
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B. Flight 


Reconstructions for the roll, pitch, and yaw axis moment coefficients for flight data are shown in Figs. 
13 through 15. As with the simulation maneuvers, both the single-axis and multi-axis maneuvers were 
performed separately. Unlike in simulation, however, the environment is not as controlled as in flight, and 
there are differences in the two sets due to expected variations in the nature of the test points. Thus for 
the flight coefficient reconstructions, the single- axis results are plotted separately from the multi- axis results 
in order to show legitimate comparisons. All reconstructions are shown for the same surface deflection 
as a representative case: 5s2 versus Sa2 for the lateral-directional coefficients, and Ss2 versus Se 2 for the 
longitudinal coefficients. 

At 6 deg Trim Alpha, Rep 1 At 6 deg Trim Alpha, Rep 1 




Figure 13. Coefficient of rolling moment for single-axis and multi-axis models, reconstructed from parameter 
estimates on flight data at 6 deg trim angle of attack. 


At 6 deg Trim Alpha, Rep 1 At 6 deg Trim Alpha, Rep 1 




Figure 14. Coefficient of pitching moment for single-axis and multi-axis models, reconstructed from parameter 
estimates on flight data at 6 deg trim angle of attack. 


The flight data reconstructions show the baseline, single-axis model performing very well, as expected. 
The multi-axis model exhibits noticeably degraded accuracy in the reconstruction, but the trends in the 
observed time histories are in fact captured quite well. Overall, it appears that the multi-axis model does 
perform close to the single- axis model in terms of moment coefficient reconstruction, achieving fit errors of 
roughly 1-5%. 

Representative plots of the reconstructions of aircraft wind angles and rotational rates are shown in Figs. 
16 through 19. 

The calculated angles of attack for both single-axis and multi-axis models, in Fig. 16, are passable 
but weak. The multi-axis model does appear to be slightly worse, but the difference is not significant. 
The calculated angles of sideslip exhibit a better comparison. In Fig. 17, the single-axis model is seen to 
very accurately reconstruct the observed angle of sideslip; the multi-axis model performs less well, but still 
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At 6 deg Trim Alpha, Rep 1 


At 6 deg Trim Alpha, Rep 1 




Figure 15. Coefficient of yawing moment for single-axis and multi-axis models, reconstructed from parameter 
estimates on flight data at 6 deg trim angle of attack. 


At 10 deg Trim Alpha, Rep 1 



(a) Single axis model, 10 deg trim angle of attack. 


At 8 deg Trim Alpha, Rep 1 



Time, sec 

(b) Multi-axis model, 8 deg trim angle of attack. 


Figure 16. Angle of attack for single-axis and multi-axis models, integrated from force and moment coefficient 
reconstructions on flight data. 


At 8 deg Trim Beta, Rep 1 



At 6 deg Trim Beta, Rep 3 



(b) Multi-axis model, 6 deg trim angle of attack. 


Figure 17. Angle of sideslip for single-axis and multi-axis models, integrated from force and moment coefficient 
reconstructions on flight data. 
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At 8 deg Trim Alpha, Rep 1 



At 10 deg Trim Alpha, Rep 3 



(b) Multi-axis model, 10 deg trim angle of attack. 


Figure 18. Pitch rate for single-axis and multi-axis models, integrated from force and moment coefficient 
reconstructions on flight data. 


At 6 deg Trim Alpha, Rep 2 



At 8 deg Trim Alpha, Rep 1 



(b) Multi-axis model, 8 deg trim angle of attack. 


Figure 19. Yaw rate for single-axis and multi-axis models, integrated from force and moment coefficient 
reconstructions on flight data. 
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manages to capture the important features of the motion. Similar conclusions may be drawn from observation 
of the pitch rate reconstructions in Fig. 18. No significant differences in the quality of reconstruction, 
however, may be readily drawn from the plots of calculated yaw rate in Fig. 19. 


V. Concluding Remarks 

Multi-axis methods provided parameter estimates that were generally similar to single-axis methods in 
terms of trends, but individual parameters were often quite different by a typical average absolute error of 
roughly 15-20%, with decreased statistical significance. Regression on flight data indicated that a multi-axis 
model exhibited an increase in overall fit error of roughly 1-5% for the linear regression estimates of the 
force and moment coefficients with respect to a single-axis model. These errors were calculated as applied 
to flight data designed for each method. Multi-axis reconstructions of aircraft state variables were shown 
to be noticeably, but not significantly, poorer than single-axis reconstructions of the same wind angles and 
rotational rates. 

Since the multi-axis maneuvers captured similar parameter estimates as both longitudinal and lateral- 
directional maneuvers combined, the number of flight-test points for the inner surfaces could have theoreti- 
cally been reduced by 50%. This reduction, however, came at the cost of moderately degraded accuracy in 
individual parameter estimates and slightly degraded accuracy in overall model fit error. Additionally, re- 
moval of half of the required test points will not necessarily reduce the overall flight time by the same amount, 
as the overhead involved in operational considerations (including getting on-condition, fueling, turns to stay 
within boundaries, et cetera) does not scale linearly and is usually highly dependent on the design of the 
test matrix. 

The inner surfaces of the X-48B Blended Wing Body, on which this study was focused, are relatively 
straightforward in terms of function and motion; application of multi- axis techniques to the more complex 
outer surface clamshell motion could result in different performance for multi-axis methods. Future parameter 
estimation research at the NASA Dryden Flight Research Center (Edwards, California) may delve into these 
problems of interest for hybrid- wing-body type aircraft. 
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Problem Statement 




• Hybrid-Wing Body aircraft such as 
the Boeing X-48B BWB are 
potentially the future of highly 
efficient air transport 

• New systems will possibly require 
new ways of approaching the 
system identification problem 

- Math is the same, but... 

— New configurations = new issues 

• Identifiability! 

- Many co-planar surfaces 

- Complex outer surface motion 
(clamshell) 

• This paper: tradeoff on accuracy 
when using asymmetric 
maneuvers to attempt to 
estimate multiple axes at once 
(time savings) 

Multi-axis versus single-axis... how to they compare? 



NASA /Tony Landis (ED06-0198-37 ) 





Surface Definitions 


• Surfaces 1: elevator 
motion 

• Surfaces 2-5: elevon 
motion 


• Surfaces 6/8 and 7/9: 
clamshell motion 

— Aileron-like for roll 

- Split motion for drag 
rudder mode 

• Winglet rudders 



Supermaneuver 



• Individual doublets are 
spliced together into a 
single data file for 
combined analysis 

• 2 sec time frame skip to 
reset integrators 

• Single-axis (pitch): 
symmetric motion 

• Single-axis (roll): 
antisymmetric motion 

• Multi-axis (pitch and roll): 
asymmetric motion 


Pair 1 

Pair 2 

Pair 3 

Pair 4 

Pair 5 






0 10 20 30 

Time (sec) 


Data Sources 




• Simulation 

- Boeing nonlinear sim v4. 3, 

VMS v4.3.1, aero model 
20100527, implemented in 
Simulink 

- PID/OBES and control 
allocator modified by DFRC 
for maneuver definition and 
forcing open-loop 

• Flight 

- Flights 65, 66, 67, 85, 89 

- Flight data filtered though 3 rd 

Order Butterworth filter, 5 Hz nasa /T ony Landis ( edio-oo56-32) 

corner frequency 




Linear Regression for PID 



• Model form is assumed 
to be linear (though 
regressors may be 
nonlinear) 

• Use statistical tools to c , = c lo 
assess model quality 
and parameter 
significance 

• Estimate uncertainty 
from Cramer-Rao 
bounds 


Outputs = Regressors * Parameters 

\ / X 

* = xe 

e = (X T X)~ 1 x T z 
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J2( z i ~ m) 2 
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Statistical Analysis 



• Analysis of coefficient of 
determination, fit-error, 
and significance suggest 
a priori that the multi- 
axis model is likely to 
perform worse than the 
single-axis model 

• Confirms intuition 
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Structure 
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Single- axis 

92.10% 

27.80% 


Multi-axis 

90.66% 

30.36% 
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31.50% 
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Selected Aerodynamic Derivs 




• Sim SA 
□ Fit MA 
A FltSA 


[] 



6 7 8 9 10 11 

Angle of attack, deg 


• Sim SA 
□ Fit MA 
A FltSA 


Derivative 

Diff. 

C Na 

16.78% 

Cmc 

18.82% 

Cyp 

82.98% 

c „ 0 

13.51% 

c w 

48.45% 
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• Sim SA 
□ Fit MA 
A FltSA 
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Roll/pitch derivs. (Sfcs 2,3) 
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• Sim SA 
□ Fit MA 
A FltSA 



* Sim SA 
□ Fit MA 
A FltSA 



— • — Sim SA 
□ Fit MA 
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• Sim SA 
□ Fit MA 
A FltSA 


Less Effective p More Effective Less Effective p More Effective 
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Roll/pitch derivs. (Sfcs 4,5) 




— •—Sim SA 
□ Fit MA 
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—•—Sim SA 
□ Fit MA 
A FltSA 



• Sim SA 
□ Fit MA 
A FltSA 



— •—Sim SA 
□ Fit MA 
A FltSA 


Angle of Sideslip Coefficient of yawing moment 



Time History Reconstructions 



At 6 deg Trim Alpha, Rep 1 



Time, sec 


At 6 deg Trim Alpha, Rep 1 



Time, sec 



N 


At 8 deg Trim Beta, Rep 1 At 6 deg Trim Beta, Rep 3 




Time, sec Time, sec 



Conclusions 

• Trends in multi-axis parameter estimates were 
quite similar to single-axis 

• Typical absolute average difference in 
individual parameters: 15-20% 

• Typical multi-axis method increase in model fit 
error: +1-5% 

• Multi-axis method exhibited greater CRB/ 
scatter at higher angles of attack 

• Multi-axis method (as studied here) could 
have reduced flight test points by 50% 




Questions? 





